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In the present paper one-dimensional two-component atomic Fermi gas is considered in long-wave 
limit as a Luttinger liquid. The mechanisms leading to instability of the non-Fermi-liquid state of a 
Luttinger liquid with two-level impurities are proposed. Since exchange scattering in ID systems is 
two-channel scattering in a certain range of parameters, several types of non- Fermi-liquid excitations 
with different quantum numbers exist in the vicinity of the Fermi level. These excitations include, 
first, charge density fluctuations in the Luttinger liquid and, second, many-particle excitations due to 
two-channel exchange interaction, which are associated with band-type as well as impurity fermion 
states. It is shown that mutual scattering of many-particle excitations of various types leads to the 
emergence of an additional Fermi-liquid singularity in the vicinity of the Fermi level. The conditions 
under which the Fermi-liquid state with a new energy scale (which is much smaller than the Kondo 
temperature) is the ground state of the system are formulated. 
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I. INTRODUCTION 

The interest in low-dimensional, in particular, one-dimensional (ID), systems has been revived in recent years in 
connection with the obtaining of degenerate quantum gases in quasi-one-dimensional magnetic and optical traps as 
well as in ID optical lattices. The previous peak of interest in ID systems about ten years ago (see references sited 
in [1, 2] and in this work) was due to the development of so-called quantum wires, viz., ID electronic systems in 
inversion layers of GaAs [3]. Degenerate atomic Fermi gases have already been obtained in traps [4-7]. It should be 
noted that Fermi gases with practically any number of particles and any intensity of the interaction between them 
can be produced in traps. This can be done by using the Feschbach resonance [7]. Thus, we can expect that various 
remarkable properties predicted theoretically for ID systems both in the framework of the low-energy Luttinger 
model [8] as well as using exactly solvable models (see, for example, [9]) can be experimentally observed in such 
systems with a quasi-one-dimensional geometry. The properties of ID metals in a fairly wide range of parameters 
can be described using the Luttinger model (see, for example, recent reviews [1, 2] and the literature cited therein). 
Accordingly, such materials are referred to as Luttinger liquids (in contrast to the Fermi Hquid in the 3D case). In 
such systems, even a weak interaction leads to a qualitative rearrangement of the excitation spectrum at low energies. 
Namely, ID metals have no well-defmed one-particle excitations. The only stable excitations in the vicinity of the 
Fermi level are collective charge and spin density fluctuations (acoustic modes). These excitations are dynamically 
independent that corresponds to complete separation of the spin and charge degrees of freedom. The interactions also 
lead to a power decay of all correlation functions over large distances and times. The response to local perturbations 
is an important problem for real systems. In this case, the behavior of ID systems also differs qualitatively from the 
situation in 3D metals. It was shown in the famous publication by Kane and Fisher [10] that the potential scattering 
of right fermions by left ones (so-called backward scattering) in a Luttinger liquid with a repulsive interaction leads to 
complete reflection of excitations from the potential at low energies. The X-ray response in a Luttinger liquid has been 
studied extensively [11-13] (including the situation with backward scattering [14]). The exchange interaction with a 
spin impurity is one of the central problems for strongly correlated systems [15]. It should be noted that if the number 
of electron channels participating in the exchange interaction exceeds double the impurity spin, the system has a 
non-Fermi-liquid (NFL) fixed point, exhibiting an anomalous behavior of heat capacity and susceptibility [16, 17]. 
The Kondo effect in Luttinger liquids was studied in [18-20]. It was shown that, as in the 3D case, the problem could 
be renormaHzed to the strong-coupHng Hmit. However, two distinguishing features are observed. First, the Kondo 
effect in a Luttinger Hquid exists both for the antiferromagnetic and ferromagnetic interactions. Second, the system 
has three fixed points, two of which correspond to the one-channel Kondo behavior and one exhibits the two-channel 
(i.e., non-Fermi- liquid) behavior (we consider the impurity spin 1/2). The conditions for a stable two-channel Kondo 
behavior relative to the exchange backward scattering were obtained in [20]. 

It is important for subsequent analysis to emphasize that, in the absence of interaction between fermions, the NFL 
state associated with a multichannel (two-channel) exchange interaction is unstable to any mechanism removing the 
degeneracy of the channels participating in exchange scattering. In particular, the instability of the non-Fermi-liquid 
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(NFL) stato in tho two-channol Kondo modol upon tho introduction of anisotropy of exchange constants in different 
channels was considered in [21]. The instabihty of the same state to potential scattering of many-particle excitations 
with different quantum numbers in quantum-dimensional structures and in metals containing impurities of d- or 
f-elements was demonstrated in [22, 23]. 

In this study, the mechanisms leading to instability of the NFL state of a Luttinger liquid with two-level (pseu- 
dospin) impurities are proposed. Since exchange scattering in ID systems is of the two-channel type in a certain 
range of parameters, several types of NFL excitations exist in the vicinity of the Fermi level: density fluctuations 
of the Luttinger liquid in the charge channel and many-particle excitations generated by the two-channel exchange 
interaction in the pseudospin channel. It will be shown below that allowance for resonant scattering of ID fermions 
(along with their potential backward scattering) from many-particle excitations generated by the two-channel ex- 
change interaction leads in our case to the emergence of additional narrow Fermi-liquid resonances in the vicinity 
of the Fermi level (even for very weak backward scattering) and to the instability of the NFL state in this sense. 
A transition from the NFL to the FL state is accompanied by the emergence of a new small energy scale and, as a 
consequence, by an anomalous increase in the density of states at low energies. 



II. IMPURITY-FREE ID ATOMIC FERMI GAS AS A LUTTINGER LIQUID 

1. Fermi gases in a trap are formed by fermion atoms of mass m with two intrinsic states [4-7]. The number of 
atoms in each state is the same. If the atoms are cooled to a temperature below the Fermi temperature Tp, they 
form a degenerate Fermi gas. The system can be treated as effectively one-dimensional if the characteristic energy of 
longitudinal motion is much smaller than the characteristic separation lo± between transverse quantization levels. In 
other words, the condition Ef <C uj± {h = 1) must be satisfied. At low temperatures, only s-type collisions are possible 
between Fermi atoms with different spins. For this reason, the interactions are characterized by a single parameter, 
scattering length a <C l±, where l± = (l/mw^)^/^. The effective ID interaction can be represented in the form of 
short-range potential with a characteristic value of <? = 2TTa/ml\ [4, 5]. Taking these limitations into account, we 
write the Hamiltonian 

H = Y.j dx^t{x)[^-dl + Ve,t{x)]^a{x)+9 J2 J da;vl/+(x)M/+(x)M/,(x)vl/,.(x) (1) 

In a box of length C with periodic boundary conditions "^aaix) = (l/vT) J^k ^^'^^'^ka (it should be recalled here that 
the relations N oo, L oo, {2'kN/L) = kp hold in the continuous limit, N being the number of fermions and 
kp the Fermi momentum), we have 

k(7 ki,k2,q 

In a ID system the Fermi "surface"is composed by two points, ikp- In the non-interacting case one has kp = nhn/2, 

n = {N-i + Ni)/£. Then, let us introduce the left(L)- (around — fc^) and right(R)- (around +kF) moving fermion 

operators aa,a{k), a = ± = L, R and the corresponding densities pa,a{(l) = 'l2^a k+q a'^a,k,<7- If we suppose, that (1) all 

fe 

the states up to the Fermi energy are occupied and (2) free dispersion relation is linear, Ska = cxhvpk, vp = 'Khn/2m, 
we obtain 

[paai-q),Pa'a'iq')] = Sa,a' Sq,q' — ; [Hq, Pa{q)] = qapa- 

Ztt 

With taking these relations into account, the kinetic energy takes the form 

27TVF 



In turn, the interaction Hamiltonian can be written as 

Hint = J2 [j dx(^gr' P<y,a{x)pa'a'{x) + Pc.,a{x)paa'{x)^ +91"' j ^(x)4'+ _^(x)*a,_„(x)*a',_<,(x) 
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For 5i > one renormalizes to the fixed line g* = 0, 53 = 32 — 5i/2 and the fixed point Hamiltonian is a Luttinger 
model [1,2]. In this case 51-interaction is irrelevant. For spin independent interactions one has (72c 7^ 52s = and 
K* — 1. These relations are generated by the equations 53, c — 52*^'^ + 91"""'^ 92. s — 92^'^ ~ 92'" ■ what follows, 
the subscript * is omitted. 

The Hamiltonian + with i/^'^j = Hint{92, 94) is equivalent to that of Luttinger model. To obtain the 
standard form of the Luttinger Hamiltonian, the densities pa,a-ix) are written by means of canonically cojugate boson 

fields (j>, and 0, Q, 0, 0: = (1/V8^)(9,0e - aH, + a{d^(j)s - aH,)). Here H, = d^d,, [(l),{x),<j),,{x')] = 0, 
[ni,(a:), n,y/ (a;')] — 0, [(j>,y{x),Il,y' {x')] — i6^^i,iS{x — x'), v — c,s. The fields 0c, s are determined by charge and spin 
fluctuations respectively: pc — dx(f>c/VT^, Ps — dx<f>s/V^- 

It should be emphasized that in an atomic Fermi gas, "charge "fluctuations correspond to fluctuations of 
the average density of the gas, while "spin"fluctuations describe the fluctuations of the relative density 
at two levels corresponding to intrinsic states of atoms. In what follows we use the standard terminology of 
a Luttinger model. 

The Hamiltonian Hq + H^^. takes the form 



,^ + 54y_/52y. I.VF+94~2g2 



(3) 



TT / V TT / V T^^F + ff4 + 2^2 

In the Luttinger model fermion fields ^^aaix) are defined as (see, for example, J3|; as well as 

^c.„{x) = exp [ - iV4^<^c.„{x)] ; $„^(a;) \a((l>c{x) + a(ps{x)) + (Oc{x) + a6',(a;))l , (4) 

rjaa are the so-called Klein factors. 

For a pure Luttinger liquid, the Green function Gaa{x; t) was obtained in [2^. The spectral functions and, accord- 
ingly, the density of states PQaa{s) were calculated almost twenty years later in [2^, ji^. In particular, the total 
spectral function (density of states) for a pure Luttinger liquid in the long-wave Hmit is defined as 

Poaa(,e) = — — - — - — — • ( — ) ; = — ; (5) 

Here, A ~ 1, £0 ^ is the cutoff energy in the Luttinger model; 77c, ?7s are anomalous dimensions in the charge and 
pseudospin channels; and parameters Kc,Ks are defined in l^l. The energy is measured from the Fermi level. 

The Hamiltonian (pj also describes long-wave excitations in one-dimensional optical lattice with the number of 
fermions per site smaller than unity (metallic state). 

The possibility of describing an atomic impurity-free Fermi gas with the help of a Luttinger model and the methods 
for experimental observation of the spin-charge separation were also considered in [2^ . 



IMPURITY HAMILTONIAN AND SCATTERING PROBLEM 

1. Let us consider the situation when a ID lattice with a number of fermion per site smaller than unity (metallic 
state) contains localized fermions with two internal degrees of freedom. The energy of a localized state is much lower 
than the Fermi level. LocaHzed fermions do not interact with one another. We assume, however, that the wave 
functions of localized and band fermions may overlap; for this reason, a locaHzed state interacts with band fermions. 
Since a two- level atom can be described by a pseudospin variable with two values of the z component, we will refer 
to this situation as a pseudospin impurity in a metallic ID lattice. We assume that the repulsion of fermions at an 
impurity site is so strong that only one fermion can occupy this site. We will consider below a ID system with periodic 
boundary conditions. In the case of an atomic Fermi gas, the pseudospin degrees of freedom in the optical lattice 
correspond to two intrinsic atomic levels of the hyperfme structure. The Hamiltonian of such a system will be written 
in the form 

H = n + Hsc; n = Hc + Hod + Hrnu Hod = J2 -^'i"'^'^; ^'^'^ =^ ^''^ + 

(6) 

^&^ = E E (rfeT'aL.«/cW.+i?.c.); i/„ = ^ (^,>+„,d, + H.c.) , 

kk' ct^a' ,a kacr 
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where He is the Hamiltonian of ID fermions, which coincides with the Hamiltonian in the Luttinger model in the 
continuous limit, H^d is the Hamiltonian of a localized level with energy Ed, rida = d'^da, Hh is the hybridization 
between band and localized fermions; Hi,s describes the potential scattering of band fermions with different quantum 
numbers (in other words, backward scattering). The term Hint from Eq. © corresponds to the repulsion between 
fermions at a deep level. 



Hint ridgUd- 



The interaction between band fermions and a localized state is generated by Hint with taking hybridization into 
account. 

The results obtained in what follows are appHcable to an atomic gas if, in addition to the relations given above, 
the inequality U < uj± holds. 

2. Excitations in a system with Hamiltonian H = H + Hsc are completely determined by the Green function of 
complex argument z, Gd{z) —< d\{z — H)~^\d >, which can be evaluated using the equations of motion. The system 
of equations of motion for Gd{z) has the form : 

g^f^-\z)Gd{z) - Vd%Gdo.{p;z)^h 

poL—LM, 

G-^\p- z)GdL{p; z)-Y, T^J^GdRip'; z) ~ V^^d^diz) = 0; 
p' 

G-i(p; z)GdR{p; z) -J^TpfGdLip'; ^) - V^i^Gdiz) = 0; 
p' 

Here Qf\z) =< d\[z - {Hod + Hint)]~'^ \d >, Ga{p;z) =< apa\{z -n)-^\apa >, Gda{k:,z) =< d\iz ~ n)-^\aka > are 
the Green functions taking into account all interactions, but disregarding scattering. In expressions 10 and below, 
the spin indices are temporarily omitted. In the case of separable matrix elements T^^" and/or those exhibiting a 
weak energy dependence the solution for Gd{z) has the form 

1 ~ hsc[z)yd(z) ^ 



(8) 



S^^^) = Y.W:\p)G^{p-z) = W^^Hpf) I de^ ^ W^\pr)E^{z), 

J Z £ 



P 

here Gdiz) is the Green function of impurity degrees of freedom disregarding scattering processes associated with 
Hbs- However, by definition, this function includes all interactions between band and localized fermions, which are 
generated by the repulsion Hint of fermions at an impurity site. The self-energy functions S°^(z) are written in the 
form of spectral expansions of many-particle Green functions for ID fermions; W^^ are the products of matrix elements, 
which are defined by the superscripts of the self-energy functions; and Pai^) is the density of states corresponding 
to the many-particle excitation spectrum taking into account the interactions. This spectrum is determined by the 
Green function Ga{p, z). The complete Green function in relations JHJ has singularities of two types. Functions Gdiz) 
contains singularities generated by the interaction between band fermions and a localized state. The denominators in 
relations l|Hl appear due to scattering of various types of many-particle excitations from one another. It follows from 
the expression for Hsdz) that the system experiences, first, potential scattering of right fermions from left ones (the 
term containing {z)!]]^ (z), but taking into account the modification of the density of states due to the exchange 
interaction and, second, the resonance scattering of ID fermions from many-particle excitations, which are described 
by the terms with the Green function Gd{z)- It will be shown below that in both cases scattering may lead to the 
emergence of additional singularities (namely, simple poles) in the Green function in the vicinity of the Fermi level. 
The poles correspond to Fermi-liquid excitations. The position of poles Zr ~ Sr + i^r is determined from the solution 
of the equation 

l-S,,(z,)-^rf(z,) = 0. (9) 
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Thus, to solve the scattering problem, we must know the many-particle low-energy density of states for ID fermions 
(disregarding scattering) as well as the Green function Gdiz) for the resonance level. 

3. In the present paper, assuming that the interactions of fermions with one another and with the impurity are 
strong, we first determine low-energy many-particle excitations, and then take into account the scattering of these 
excitations from one another. In other words, first the problem with Hamiltonian H is solved and Gd{z) is determined 
with the help of the methods used in the problem with J7 ^ oo. 

Because the energy U is dominant, it is essential to properly treat correlation on site, and the method of the 
equations of motion correctly accounts for these on-site correlations This method gives the following expression 
for the Green function of band fermions 

Gaaa{k, k' ; z) = 5{k ~ k')gaa{k\ z) + Gaaik] z)Taa{k, k' ] z)gaa{k'; z), (10) 

where Gaik; z) is the Green function of a pure Luttinger liquid, and the scattering matrix has the form 

TUk.k'-z) = V^^GaizW^,. g,iz) = g^r'\z) - ^S]o^J(z), 

a 

Supposing that in the limit of large values of U the main effect of this interaction is the many-particle resonance 
at the Fermi level, we use the following relation for low energies 

Qa\^) = E So'7(^)'«< (He + H,, + H,, + H^^jw >, (11) 

a 

here Sj^(z), H^J are defined with poaais) corresponding to a pure Luttinger liquid, H^x is the exchange scattering 
of the ID fermions by two- level (psevdospin) impurity. In accordance with the relation ljll|l . to calculate Qd{z) in 
what follows, we represent the exchange interaction in the form of the model of a resonance level. 
By means of Eq. ifTTijl , the density of states is defined as 

PQcr(£) = T7r~"^Im 5pGaaCT(fc, A:'; e) ^ poaa{e) ^ ■n^Hm SpQd{e)'^\V^^j\^Ql^„{k;£) = poaa + Pdaa, (12) 

k 

where the minus and plus signs correspond to £ > and £ < 0, respectively. The relation itTTll impHes that only the 
density of states pdaa generated by the exchange scattering enters the self-energy Ssc(-z) in Eq. ||H||. 



THE RESONANCE-LEVEL MODEL FOR THE EXCHANGE INTERACTION 

In one-dimensional systems the exchange interaction takes the form H^x — Hex ^ + Hex ^ , where 

HiP-T. E E-^'*-(o)^^v*-'(o)s^ i?if)=E E J2JBHA0K.'^a'a'm m 

a i—x,y^z a^cr' a^a' i—x^y.z cr,o"' 

Here, S is the impurity pseudospin, and are the Pauli matrices. Hamiltonian Hex corresponds to the back-ward 

( p\ 

exchange scattering emerging in the presence of anisotropy of scattering channels. Hex describes the two-channel 
exchange scattering. 

The interaction in the charge channel of a Luttinger Hquid is characterized by parameter (see expression 
(??) below). It was proved in [20] that, in the case of repulsive interactions {Kc < 1), the small anisotropy (J_b) 
of the exchange scattering channels in the vicinity of the fixed point of strong interaction in J is relevant only for 
1/2 < iiTc < 1. In this case, the results obtained in [19] are valid (it was found in [19] that the two-channel Kondo model 
is absolutely unstable). However, for Kc < 1/2, the two-channel Kondo behavior is stable to the exchange anisotropy 
of the channels. In this case, however, one more mechanism violating the channel degeneracy (weak resonant and 
potential scattering of many-particle excitations) in the two-channel exchange interaction was discarded. We will 
consider this mechanism here. Using the fact that the exchange anisotropy of the channels is irrelevant for Kc < 1/2, 
we consider the situation when the impurity pseudospin has a symmetric coupling with adjoining lattice sites so that 
Jb = 0. It will be shown below that Hamiltonian itT^ll in this case can be reduced to the model of the resonance level 
(this model was constructed in [3| for the two-channel exchange scattering in metals disregarding an interaction 
between band electrons). 
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For ID system, in order to represent the two-channel exchange interaction in the form of the model of a resonance 
level, we introduce the boson representation Q of fermion fields "^aaix). Hamiltonians H^, in Eq. (3) lead to the 
canonical (for bosons) form 



with the help of the following redefinition of boson fields: 4>i, — \/Ky ■ (j)'^; li^, = [l/^JK^) ■ IV^. In this case, phases 
<i>as in relation Q are transformed to 



/Kc VK, 

Substituting the fields (j)^, 9^ into Eq. ltT3|l and taking into account relation itTljl . we obtain 



V Y J\. c ■'J X — 



ttK, 



■{dxO'g)x=o- 



(14) 



Thus, in the general case, the exchange interaction in a Luttinger liquid has the form 



cos (V4^-0'.(O))e-V(4-/if.)-eUo)+^.c.] + j^s^,J^^Q^e%^^ 



(15) 



It should be emphasized that the expression for the exchange interaction contains only pseudospin fields. To 
reduce expression II 1511 to the Hamiltonian in the model of a resonant level, we carry out the following transformations. 



1. We introduce fields $l,_r, instead of 9'^ 



^^±^- ^' = ($^ + $«) . 



2. For convenience, we replace the right field $_r by the left field <i>^ : ^'l{x) ~ — $fl(— .t). 

3. We introduce the symmetric and antisymmetric fields ^s,A '■ = — — ; $r = 



V2 



V2 



As a result of these transformations, Hamiltonian lfT5|l assumes the form 

Ji 



2'Ka 



27r 



(16) 



Field $s from the transverse part of the interaction is eliminated by rotation about the axis: U = e2'*sS". For 
gi = (or, which is the same, for Kg = 1), the kinetic energy in the pseudospin channel can be reduced to the form 
corresponding to free boson fields so that 



H, = 



VF 



(17) 



Since UHgU ^ = Hg — A{vf lT^){dx'^ s)S'^ , we obtain, using the transformation U{Hs + Hex)U ^ , the Hamiltonian 

is'] 

Hs = Hs + Hex in the pseudospin channel in the model of a resonant level; here. 



H. 



Ji 



V-KO) + ^^(0)1 (d+ - d)) + -^(J, - %vf) ■ ^J(0)^s(0)(d+d - 1/2) = Ha + Hs 



2'Ka 

g-2i$s,A(a:) 



27r 



'2'Ka 



; S+=d+f]] S*^ = (d+d- 1/2); f = I. 



(18) 



It should be noted that phases <^s.a{x) describe the symmetric and antisymmetric fluctuations of pseudospin density. 
The total Hamiltonian of a Luttinger liquid with a two-level impurity has the form 



H^He + Hs + Hit} + Hsc ^Ho + M$ + H,. 



(19) 
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If we disregard Hsc, Hamiltonian H19|l gives excitations of the Luttinger liquid with velocity Uc in the charge channel. In 
the pseudospin channel, excitations i^s a determined by the two-channel exchange interaction, generating, among 
other things, a resonant many-particle level for low energies. A remarcable feature of the resonance-level model ifTsjl 
is that the hybridization and interaction are performed via different channels. This is its essential difference from the 
resonance-level model for the single-channel Kondo scattering, in which both of them are in the same channel. 

2. The model Ijl8|l can be renormalized to the limit of strong coupling with the fixed point being on the line 
Jz = {Jz — 8vp) — [21] (the Emery-Kivelson line). On this line, the impurity degrees of freedom are hybridized only 
with the field associated with antisymmetric pseudospin fluctuations. For T — 0, the Green function in the vicinity 
of the Fermi level has the form |2JI 



To - To 



e + iTKsigne e + i^signe 



(20) 



In the strong coupling limit, Tr- ~ Tk, Tk being the Kondo temperature. The second term in relation Ij2()|l emerges 
due to the fact that half the impurity degrees of freedom in the two-channel Kondo model are not hybridized with 
collective variables of band fermions. Since Hamiltonian ifTsll does not conserve the number of fermions, Gd{z) has 
nonzero anomalous matrix elements dd > and >. 

For the problem discussed in this paper, the form of Gd{z) beyond the Emery-Kivelson line or, what is the same, 
for a finite interaction in the S-channel, i.e., for Jz — Svp 0, is important. To obtain a solution for this case, let us 
use the technique that was previously applied to the well-known problem of x-ray absorption in metals [29]. First, the 
Hamiltonian Hqs + Hs, here H^s — ivp ipg{x)dxilJs{x), of Eq. ifTsjl is diagonaHzed. To this end, we introduce 
boson operators bsk — k^^^^Ps{k), b'^f, = k^^^^ ps{—k), where Ps{k) are density operators 

^ ku—k ^ ko 

^«(^) = ^t{<i)^s[q + k) : ps{-k) = ^ ^tiq)Mq ~k), k> 0, 

13=0 q=k 

ips{k) are Fourier components of fields ips{x), and the cut-off takes place at ko ~ a~^. Using operators bsk, ^^j., we 
write the Hamiltonian Hqs + Hs as 

Hos + Hs = vfJ2 ^Kk^sk + K{d+d -\)Y. ^'^'(^i + ^-fe)' (21) 

fc>0 fc>0 

Here \z = {1/2tt){Jz — S£f)N~^/'^ . This Hamiltonian is diagonalized to become vp X]fc>o ^^tk^sk by the canonical 
transformation 

U = e^plxzPo{d+d-^)Y,k-'/Hbsk~btM , po-^SpK 

\ k>Q ) 

In this operation, the Hamiltonian Ha is transformed to 



Ha = + ^Aimi^ - d) + A{d+d- i), 

d+ = Ud+U-^ = exp Ixzpoa k-'/^{bsk - 6+fe) J d+ = Uod+ , 

\ k>0 ) 



(22) 



A = — ej — XIpo is the "polaron shift"due to the screening interaction in the S-channel. With due account of Eq. 
12211 . the Green function Gdit) of resonance level is 



Uolt) is derived from Uo{0) using the substitution bsk 6sfce'^'=*. Here (...)_d denotes averaging over the states of the 
diagonalized Hamiltonian Hqs + Hs, G^iit) is the Fourier transform of the function given by Eq. Ij2()|l . The averaging 
is performed in the convential manner using the relations 
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e^e^ = gA+B+(i/2)[A,B] 



,(l/2)(F=(6+,&)) 



where F is an arbitrary linear combination of boson operators. As a result, we find that, at large times ept ^ 1, the 
function in Eq.lj^Sl 



Consequently, we obtain in the energy representation the expression 



To 



A + iT 



K 



£ - A + iVK 



£- A 

So 



£- A >0, 



(24) 



Here, A ^ e^^, T{x) is the gamma function, as = {5s /tt)"^ vi 5s ^ is the phase shift due to the screening interaction 
in the S-channel, A = — £j. The cut-off parameter £o '--^ £_f since the velocity of excitations in the collective channels 
is vf, as follows from the Hamiltonian Ijl8|l . 

It should be noticed that the power-law dependence 1I2H1 occurs when the condition Tk S> \ej\ is failed. In this 
condition the Kondo temperature is defined on the Emery-Kivelson Hne, Tk ~ £f exp[— 1/( J_Lpo)]- Using the expres- 
sions for Tk, \sj\, we find that the power-law dependence 1I2H1 takes place under condition {J±/eF) *C 1/ ln(£F/| Jz|), 
Jz = (l/27r)(J2 — 8wf). When the inverse inequality is satisfied, one may use the expression Il2()ll on the Emery-Kivelson 
line. 



CALCULATION OF THE DENSITY OF STATES AND SELF-ENERGY FUNCTIONS 

1. The Green functions of impurity degrees of freedom for low energies are defined by relations Ij2()|l or Ij24|l . Thus, 
to determine the many-particle density of states for low energies in the presence of a pseudospin impurity, one needs 
to substitute these Green functions into ltT2ll . 

First, we must know the quantity X;fe l^kdl'^GlAkie) ~ \VC^d\^T,kSL(k;e), where X;fc e) = ~di:oaa{e)/de. 

Functions T,oaa{£) are defined by expressions © with the density of states paa{£) — Poaa{£) from @; namely, 

So..(-) = j de ^-^ = ^^^-^[j de + Jde -^J = A„( Jo- + Jo+), (25) 

-eo 

Here, z± = e ± ij, e > 0. For £ < HMeen T,g~l{z) = — Eq^^(z). Integrals Jo^(z), which are Gilbert transforms of 
the density of states for a Luttinger liquid, are defined as 

»7c + l 

iVc + 1)z z± 

where F(a, l3]j; x) is a hypergeometric function. Since we consider here jmly the case when £o/|-z| ^ 1, we can use 
the transformation formulas for a hypergeometric function with |a;| ^ 1 ,30] and obtain the following expressions: 

r(/?)r(7 - a) r(a)r(7 - p) 

We have used the familiar relations r(l) = 1; r(r/c + 1) = r]cT{ric); r(l — ric)T{ric) ~ it/ sin(7r?7c). Thus, in long-wave 
limit, for 7^0, we obtain 

vW _ / 2z7rsin[(7r/2)?7c] exp[i{n/2)rjc] \ „^ , . 

^ sm(7r?7c)r(l -I- 77c)eo ^ 

Since F{1, 1; 2; ^x) = ^x~^ ln(l ± x), expression (j2(ill for ijc = corresponds to the principal term in the expansion 
in parameter x = (£0/^) ^ 1- Using the expressions 

T.SL{k;e)^y-^§^, £>0; T.SLik;e) = -^^^^-^^^, £ < 0, 
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we obtain with the help of formulas (j25|l and (j2(ill 



— 2« Gxp[i(7r/2)77c] \ / TTr]c sin[(7r/2)77, 

r(i + ^c)4+"^ 



(27) 



For the Green function Ij2()ll (or, what is the same, on the Emery-Kivelson line), substituting the expression Ij27|l 
into Ijl2|l . we arrive at the following expression for the impurity contribution to the density of states: 



2V^ricSm[{TT/2)'nc 



r(l + ryc)£j^''" sin(7r77c) 



cos[(7r/2)?7c 



sin[(7r/2)77c 



^ K 



K 



+ cos[(7r/2)r7j • \e\ 



(28) 



Here, we have introduced the notation V'^ = IVJ^^P- For low energies, the most singular term in the density 
of states is Pda[\£\) ~ |e|'''=~^. Pay attention to the fact that this contribution is due to the term of the form 
R.e5d(£^) • ImX^fc ^aCT(^;£) in R-efJd(e) we take the part corresponding to nonhybridized degrees of freedom (second 
term in relation Il2()|l l. While deriving formula Ij28|l . we took into account the fact that the term with the i5-function in 
Init?rf(e) makes zero contribution to the density of states for rjc ^ 0. Thus, we see that, in contrast to noninteracting 
Fermi gas, excitations in a Luttinger liquid considerably modify the peaks at the Fermi level, which are generated by 
the interaction with a two-level impurity. The singularity in the impurity density of states is enhanced for 77c < 1; on 
the contrary, it is suppressed for rjc > 2. 

2. The contribution to ^a{z) from the impurity term pda(e) to the density of states is determined by the expression 



+00 



E.(z)=/ de^^ /cie^^^I]W(z+) + I](-)(z_). 



z — e J 2+ — e J + e 





(29) 



Taking into account expression (I28|l for the impurity density of states, we see that Sq(z) contains the integrals of 
three types, which are defined in the complex plane: 



r(T) 



de 



^ At) 



-j-CXD 



; n 



de 



£ At) 



; 4 



de 



Evaluating integrals /2^'' (which are singular for small z) and substituting the obtained expressions into (I29|l . 
we obtain 



(1) for Pda 



,(-1) for " ±"; Ai 



2m]c coa{nrjc/2) sm{'Ki]c/2)V'^ 
eg""*"^ sin(7r77c)r(ryc + 1) 



(30) 



(2) for Pdc-^ler^-'; sW=A2-z^=-^-P±; 
P± = -e-"-, (+1) for " ± "; A,- ^-VcSin'i-Vc/2)V^ 



Tk sin^(^?7e)r(rye + l)e^ 



(31) 



Integrals /g^^ (z) are defined by the same formulas that were derived while determining Sq^^ . 

3. Beyond the Emery-Kivelson line, expression 1I2H1 defines the retarded Green function Q^{e), e = e — A. Let us 



consider the case when 



J2k^aai^T^) J cis a function of energy varies much more slowly than Q^{e). The 



definitions of these two quantities imply that this is possible for ijc ^ 013. In this case, we can set '^f. Q^^{k; e) 

For the nonhybridized degrees of freedom (second term in the expression it can easily be verified that the 

density of states for has the form 



(32) 
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The retarded self-energy functions are defined as 

+ CXD OC A CXD 

(.) = / deJ^ = fde^ + J ds^ = fde\^ + , (33) 

J z+-e J z+ -e J z+~£ J - - - I- 

— oo ^ — oo 



z+ — £ z+ + e 



where z^ = z+ — A; Rez+ = £ > 0. Using formulas (I33|l and the expressions for the integrals derived above, we 
obtain 

Sf(z)--Q4— -(^ ; Qd = QdCig{T:as). (34) 



z\ 



Since as < 1, the above contribution of the impurity density of states to (z) is also singular for small 

FERMI-LIQUID RESONANCES IN THE VICINITY OF THE FERMI LEVEL 

Let us now consider the Fermi-liquid singularities in the total Green function which are associated with the 
scattering of many-particle excitations, described by Hamiltonian Hsc in expression (EJ. In the presence of resonance 
and potential scattering, the position of the poles is determined from the solution of Eq. lO- Substituting Sgc into 
Eq. lO and assuming that all matrix elements appearing in are determined by their values for k = kp, we obtain 
the equation 

1 - T^^L{Zr)^R{Zr) ~~ T/^T^d (2r)SL (z.)I]fl (z.) [2 + r(EL(z,) + SflX^r))! = 0, (35) 



where V and T are the matrix elements of the resonance and potential scattering and Sq are the self-energy functions 
calculated in (l3?11l , ll3Tll , and l(3ijl . The three terms on the left-hand side of Eq. ll35)l , containing self-energy functions 
Sa, describe different scattering processes of many-particle excitations. The term proportional to corresponds to 
scattering of right fermions from left fermions taking into account the fact that the density of states contains the 
impurity contribution pda- The terms proportional to and T^, describe two possible processes of resonant 
scattering, i.e., scattering of charge and pseudospin densities from many-particle impurity degrees of freedom, which 
are described by the Green functions Gd{zr) defined in ll2(11l . ll2Hl . 

1. On Emery-Kivelson line, the impurity density of states is defined by expression ll28ll : accordingly, the poles 
in the vicinity of the Fermi level are generated by the contributions to Eq,, which are most singular in a certain 
range of parameters for small |z| ^ Tk- We will first consider the solutions to Eq. Ij35|l . which are associated with 
scattering of nonhybridized excitations with band collective excitations of impurity degrees of freedom. In accordance 
with expression (l2?11l . nonhybridized degrees of freedom are described by the Green function Gd{zr) = ^/zr- We begin 
with the case when the main contribution to Ea(z) comes from the term ~ |£|'''^~^ in the density of states. Retaining 
in Eq. II35|I the principal terms with small values of \z\ , we obtain the following equation for the poles: 

1 - V^T^gd{Zr)^L{Zr)^BiZr) (Si(z,) + I]il(z,.)) = 0, (36) 

For Zr±, we substitute T,^{zr±), respectively, into Eq. l(3fijl . It should also be noted that T,L{zr) = Ei?(zr) since 
= Pi?(|e|)- Let us consider solutions with Zr-. We write Zr- in the form Zr- = \zr-\exp{i(p), physical solutions 
corresponding to values ofO < f < 7r/2. Solving the imaginary and real parts of Eq. l(36j) . we obtain the following 
solution with \zr"\ ^ Tk'- 

£0 Veq/ V£o/ 7 - 377c 

where Ai is a factor on the order of unity. This solution exists when the following conditions are satisfied: 



The first system of inequalities emerges from the requirement of vanishing of the imaginary part of Eq. Il3(i|l . When 
conditions Il38|l are violated, Eq. I|36ll may have solutions due to two other contributions to the density of states. For 
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Pda ~ l^l'''' ^ and Qd{zr) = '^/zr , there exists a solution Zr- = Sr- < 0, corresponding to a localized level below the 
Fermi energy. The position of this level is determined by the energy 

Finally, for pda ~ 1^1'''= and Gd(zr) = l/^r, we have a resonance with 



which exists when the following conditions are satisfied: 

1 fT\^/i/V\ /rK\(5-2'7c)/8 

'?=<^; - • - « — • (41) 



6 Veo^ ^£0^ ^ Eo 

It can be seen that all the resonance obtained above are formed as a result of scattering of nonhybridized impurity 
degrees of freedom by many-particle excitations, which form the density of states at low energies. The scattering of 
many-particle excitations with Qd(zr) — ^l(zr± + ?rxsigner.±), which is described by the term in S^c proportional to 
T^y^, does not lead to the formation of resonances with \zr\ <^ Tk- In other words, Eq. has no physical solutions 
in this case. 

Beyond the Emery-Kivelson Hue expressions for Qd and Sa from (I2H1 and lIMjl should be substituted into Eq. 
At Re(zr) > the position and width of new resonance are determined by the formulas (we consider only contribution 
of the nohybridized impurity degrees of freedom, which is described by the second term in Eq. Ij24|l l 



/T\l/2(l-as) /\/x2/(l-as) / A N ('?o-l)/4(l-as) TT 

£0 ^£a^ ^£0^ ^£0'' 4 



here ^3 ~ 1. Since A <C £0 and rjc^ as, the resonance for rjc > I can be quite narrow with a width <C Xz/£f- 

2. In the range of parameters where Eq. (j^ffill has no solutions, the poles can be due to other scattering processes 
apart from those making contribution to S^c, proportional T^V^. Assuming that this contribution is small in the 
regions where Eq. Ij36|l has no solutions, we consider the equation 

1 - 2V^Tgd{Zr)^L{Zr)^R{Zr) = (43) 

On Emery-Kivelson line, when the main contribution to Sq, comes from the term proportional to ~ |e|'''=~^ in the 
density of states, the Fermi-liquid resonance is formed due to scattering of many-particle excitations by the resonant 
level with 

1 i 

Gd{zr) = 7 — — : r « Tt^, for Re(zr)^0. 

[Zr± + iTKSigner±j Tk 



The position and width of the Fermi- liquid resonance are determined by the solution to Eq. Ij43|l , 

/T\ l/2(2-7)c) /V^N6/2(2-7,e) / £o N 1/2(2-,,,) TT 



So \£oJ Veo/ Vr^^ ' 4(2-7?,) 

where A is a factor on the order of unity. This solution exists when the following conditions are satisfied: 

3 /r\l/6/y\ /rK\(5-2r,o)/6 



3 /r\l/6/yN /rK\(5-2r,o)/6 



The same scattering process generates the Fermi-liquid resonance with Zr- in the case when the main contribution 
to Eq. comes from the term proportional to ~ |e|'''=~^ in the density of states. In the present case, the new resonance 
exists for 



r?c<l; (- ■(- «(— • 46 



Finally, new singularities may appear due to potential scattering of right fermions from left ones taking into account 
the impurity density of states Ij28|l . The positions of the poles in this case is determined by the solutions to the equation 



12 



1 — T^Si(zr)Sij(zr) = 0. It can easily be verified that a solution exists in the case when the main contribution to Sq, 
comes from the term proportional to ~ |e|'''=~^ and has the form of a localized level below the Fermi energy. 

It can also be proved that potential scattering beyond the Emery-Kivelson Hne generates a locaHzed level above 
the Fermi energy. 

If the density of states at the Fermi level is determined by expression poa for a pure Luttinger liquid, the potential 
scattering of right fermions from left ones does not lead to the formation of additional Fermi-liquid resonances or 
levels at low energies. This case was treated in [10]. 

Concluding the section, let us prove that the inequalities for parameter rjc derived above, which determine the 
ranges for new Fermi-liquid resonances, are in accordance with the inequality Kc < 1/2. It should be recalled that 
this is the condition of appHcability of the model with Jb — 0. In particular, the conditions for rjc in ll38|l and 
correspond to the inequalities 0.13 < Kc < 0.2 and 0.11 < Kc < 0.2, respectively. Finally, the condition ijc < ^ in 
Il4(ill corresponds to values of Kc < 0.2. In this case, the mechanism considered in this study ensures instability in the 
region 1 > 77c > 1/8, which is also in agreement with the conditions derived above. 

CONCLUDING REMARKS. THE BEHAVIOR OF PHYSICAL QUANTITIES 

Several types of NFL excitations with different quantum numbers characterize a Luttinger liquid with two-level 
impurities at low energies. The results describe above indicate that the scattering of many-particle excitations of 
various types from one another leads to the emergence of an additional Fermi-liquid resonances in the vinicity of the 
Fermi level in a ID system. The conditions are dtermined, under which the fermi-liquid states with a new energy 
scale muchsmaller than the Kondo temperature is the ground state of the system. 

The type of many-particle excitations and, accordingly, the type of the phase state are determine by the following 
parameters: the value of anomalous dimensionality r]c in the charge channel of a Luttinger Hquid (or, what is the same, 
the value of parameter Kc) as well as the initial energy Ed of a deep level. This energy appears in exchange constants. 
Reducing (increasing) the depth of the impurity level leads to increasing (reducing) the Kondo temperature. Let us 
define the phase states of the system in dependence on the values of these parameters. 

In accordance of conditions Il38|l . Ij41|l . Ij45|l . and Il46|l . new Fermi-liquid resonances are absent when, at any 
rate, 77c is greater than 5/3 and/or the Kondo temperature on the Emery-Kivelson Hne is found to be very low. In the 
latter case, the impurity level is very deep and conditions of the type |zr| ^ Tk are violated and we obtain the phase 
state of the system , in which excitation of the Luttinger liquid take place in th charge channel and the excitations 
generated by two-channel exchange interaction occur in the pseudospin channel. 

The low-temperature behavior of the heat capacity C(T) is determined by the product 7T, where 7 ^ ln{TK/T) in 
the second order of perturbation theory in Jz/J± (2J|i 1^- The logarithmic dependence of C/T occurs instead of the 
constant value of 7c ^ uj^ in a pure Luttinger Hquid. The homogeneous static susceptibility behaves analogously on 
the Emery-Kivelson Hne. According to the conditions described in conclusion of section IV, the logarithmic dependence 
is transformed into the power dependence 7 ~ j^-i+as occuring beyond the Emery-Kivelson line at the crossover 
temperature Tc ~ \ej\. If the initial values of the parameters is such that the condition Tk ^ \sj\ holds, the crossover 
takes place at deepening the impurity level. 

In the phase state considered here, anomalous correlations exist at the impurity site. The divergence of the correlator 
< S~^S~ > {uj — 0,T) ^ ln(l/T) corresponds to free rotation of th impurity pseudospin. We recaH here that in an 
atomic Fermi gas pseudospin fluctuations describe the fluctuations of the relative density at two levels corresponding 
to intrinsic states of atoms. The anomalous behavior of the one-site correlator < S'^S~ > corresponds to anomalous 
increase (upon cooling) of the correlations between the occupancies for two intrinsic states of an impurity site. This 
made indicate a tendency to the formation of the superfluid stae in the relative density of two component of the 
impurity subsystem in a quasi-one-dimensional system. 

On the Emery-Kivelson line Fermi- liquid resonances exist for relatively small values of tJc and high values of Tk 
in accordance with the formulas derived above. If the value of rjc satisfiesany condition of the existence of resonances, 
but the corresponding condition of the type jz^l ^ Tk is violated, we can obtain an additional Fermi- liquid resonance 
by reducing the depth of the impurity level and, hence, by increasing the Kondo temperature. 

Expression (I42II implies that Fermi-liquid resonances exist beyond the Emery-Kivelson line for all admissible 
values of as (i.e., impurity level depths), but only for strong interactions in the charge channels of a Luttinger liquid. 

Both on and beyond the Emery-Kivelson line, transitions from an NFL state to an FL state can take place. The 
characteristic crossover temperatures are Tr jr, where are the widths of Fermi- liquid resonances. 

It should be emphasized that the Fermi- liquid resonances lead to the existence of new energy scales 7^, which are 
much smaller than the Kondo temperature. 
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